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■ Abstract. The gauge invariant electromagnetic Wigner equation is taken as the basis for 

, a fluid-like system describing quantum plasmas, derived from the moments of the gauge 

^vq | invariant Wigner function. The use of the standard, gauge dependent Wigner function is shown 

to produce inconsistencies, if a direct correspondence principle is applied. The propagation 
' (— i ' , of linear transverse waves is considered and shown to be in agreement with the kinetic 

<**> |' theory in the long wavelength approximation, provided an adequate closure is chosen for the 

macroscopic equations. A general recipe to solve the closure problem is suggested. 
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1. Introduction 

The Wigner function is the quantum equivalent of the classical particle distribution function 
and can be used to calculate average values of physical observables HI. In most cases, the time 
■^j- ■ evolution of the Wigner function is evaluated considering only scalar potentials, hence without 
the inclusion of magnetic fields. One reason for this is the considerable analytic complexity of 



the electromagnetic Wigner equation. Indeed even the electrostatic Wigner equation already 
is a cumbersome integro-differential equation which hardly can be examined except in the 
linear limit. However, the emergence of new areas like spintronics where magnetic effects 
are crucial makes it desirable to have quantum kinetic models allowing for nonzero vector 
potentials. In this situation the gauge invariance of the Wigner function should be assured 
from the very beginning in order to avoid inconsistencies, a point somewhat neglected in 
previous studies. It is the purpose of this work to stress the relevance and properties of 
the gauge invariant Wigner function (GIWF) |f3l @] |5l in connection with quantum plasmas 
problems. In addition we provide a macroscopic (moments) formulation starting from the 
electromagnetic Wigner-Maxwell system, substantially generalizing the recently introduced 
moments system derived from the Wigner- Poisson equations [6J. The resulting macroscopic 
equations are a step toward the inclusion of spin dependent variables, postponed to future 
considerations. 

The advantages in macroscopic formulations are in their relative simplicity, so that the 
nonlinear regimes are not necessarily unaccessible apart from numerical simulations. Notice 
however that our fluid approach does not imply any fluid approximations, in the sense that 
we are not supposing a large collision rate or a short mean free path for instance. If we are 
interested only in basic quantities like particle, current or energy densities, nothing forbids 
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to compute moments of the Wigner function in order to derive fluid-like equations for the 
time-evolution of these variables. The roots of the moments descriptions in plasma theory 
can be traced back to Grad [7J. The price of replacing the more detailed kinetic models by 
macroscopic models is the loss of information on kinetic phenomena like Landau damping, 
the plasma echo and many others. 

This work is organized as follows. In Section II we briefly review the definition 
and properties of the GIWF. Section III develop the corresponding fluid moment hierarchy 
equations. Section IV consider the propagation of transverse waves and the closure problem 
in this case. Section V is dedicated to the conclusions. In addition, it is included the Appendix 
A where the closure of the fluid-like system is discussed. 



2. Basic properties of the gauge invariant Wigner function 

A sensible definition of gauge invariant one-particle Wigner function / = f(r,\,t) was 
introduced by Stratonovich [3], and rediscovered by Irving [8J. Since in this work we are 
not concerned with relativistic phenomena, we write it in a non-co variant form, 

/(r ' V ' = (^h) 3 I dSQXP [l (^ + ^£ / / 2 2 ^A(r + TS,0) 

x ¥ * (r+^Ayfr-^t) , (1) 



2 / V 2 

where r and v are the position and velocity vectors and t the time. The wave function is 
assumed to be normalized to unity. In addition, h is Planck's constant divided by 2n, A(r, t) 
is the vector potential, m and q are the mass and charge of a particle in a pure state described 
by a wave function i//(r, t) . The properties to be discussed in this Section hold equally well in 
the case of mixed states. In contrast to the original definition of Wigner function [Q]| via the 
canonical momentum, the object / in Eq. (0Q) is written in terms of the kinetic momentum mv. 
The extra integral in Eq. (OQ) containing the vector potential compensates for the change in the 
wave function in a local gauge transformation. The use of a non-covariant, one-time pseudo- 
distribution renders the interpretation issues of / less obscure than in a four-dimensional 
space-time version, as stressed in Ref. |[9l- 

Naturally there are other ways to obtain GIWFs, e.g. through certain path integrals 
involving the vector potential ifTOll . However, the phase factor in Eq. (Q]) can be justified J51 
in terms of the minimal coupling principle. Moreover, as discussed in more detail elsewhere, 
the function of the phase factor is to convert any gauge into the axial gauge [5J. For our 
purposes, the choice of the form (OQ) is due to convenience as it provides a non-ambiguous 
way to calculate averaged quantities. If instead one take a GIWF in terms of a line integral 
A(s,f) • ds one introduce the further difficulty of the choice of integration path from ri to 
r 2 (cf. Eq. (2.157) of Ref. IfTOl ). 

The properties of the GIWF have been detailed in |@][5]|. Nevertheless for completeness 
we discuss some of them once again. From / we can compute the very basic zeroth and first 
order moments 

"dyf =|v| 2 , (2) 
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dwf = ^(yVV* ~ V*Vy) ~ -I VTA, (3) 
2m m 

with the interpretation of particle and current densities respectively. By construction these 

quantities are invariant under the local gauge transformation 

A^A + VA, yr^ y/exp ( l \- ] , (4) 



h 

where A = A(r, t) is an arbitrary differentiable function. 

If the starting point is the usual (gauge dependent) Wigner function 

' c >-">" = (w/" s ^ (rr) v ( r 4 «) * ( r - \ <) • (5) 

which is written in terms of the canonical momentum p = m\ — qA, one obtains gauge 
independent results for the zeroth, first and second order moments, but gauge dependent 
quantities when considering higher order moments. Of course implicitly we assume that all 
physical objects should be gauge independent. Serious discrepancies occurs when calculating 
the evolution equation for the second order moment of the usual Wigner function and the 
GIWF, as will be shown in the next Section. In all cases it is safer to work with / as given in 
Eq. ©• 

The time-evolution of the GIWF was considered already by Stratonovich but a 
particularly illuminating form to express it was provided by Serimaa et al. JH according 
to 



|| + (v + A v).l. + i[E + (v + Av)xB].A}/(r,v,0=0. 



(6) 



Here, we introduced the operators 



A _ ihq d f 1 / 2 , _ / ihx d \ 

A\ = ^- — x dzzB r+ T-,f , (V) 

m L d\ J -i/2 V m d\ J 

■V2 / ihx d \ 

M r+ (8) 



1/2 V m d\' 

, 1/2 / ihx d , 
B = / dxB r+ — ,H , (9) 



'-1/2 \ m (9v 
where B = B(r, and E = E(r, are the magnetic and electric fields respectively. The kinetic 
equation © follows from the Schrodinger equation for the wave function or, alternatively, 
from the von Neumann equation solved by the density matrix. 

As apparent from Eq. ©, the kinetic equation satisfied by / is formulated in terms of the 
physical fields, unlike the equation solved by f GD which is written in terms of the scalar and 
vector potentials ifTTI and which can be shown to be not gauge invariant, a serious drawback. 
Moreover Eq. © is almost in the form of a Vlasov equation, with two differences: the 
electromagnetic fields are replaced by E and B defined in Eqs. dH])-©; the velocity vector is 
displaced by the intrinsically quantum mechanical perturbation Av defined in Eq. ©. Notice 
that this perturbation Av vanishes in the electrostatic case. In calculating Eqs. CO)-©, it is 
assumed that the electromagnetic fields are analytic, so that the integrals are evaluated after 
Taylor expanding and then replacing r by the indicated argument r + ih ( x/m) d/d\. A further 
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difference in comparison to the Vlasov equation is that not any function / on phase space can 
be taken as a Wigner function. Too spiky functions violating the uncertainty principle should 
be ruled out. And, of course, the Wigner function is not strictly a probability distribution, 
since in general it is negative in certain regions of phase space. 

To sum up the pseudo-distribution in Eq. (Q~|) provides a practical and non-ambiguous 
recipe for a GIWF and Eq. © is the associated kinetic equation. In the next Section we 
derive a system of partial differential equations satisfied by macroscopic quantities obtained 
taking moments of the GIWF. 

3. Fluid moments hierarchy 

In spite of the apparent simplicity, actually Eq. © becomes quite complicated after 
developing the operators E and B. In practice, nonlinear problems are unaccessible in this 
formulation, specially remembering that the electromagnetic field should be self-consistently 
determined through Maxwell equations. Hence apart from linear problems this Wigner- 
Maxwell system can be helpful only by means of numerical simulations, which are themselves 
not evident due to the complexity of the system. This motivate the creation of alternative 
models capturing the essentials of the quantum plasma dynamics. 

In this context, recently [|6l a fluid moments hierarchy was derived from the electrostatic 
Wigner equation. As usual in moments theories [7], a set of macroscopic variables (particle 
density, current etc.) were defined in terms of integrals of the Wigner function. The time- 
evolution of these quantities was then deduced from the Wigner equation. No assumptions 
were made on the particular local equilibrium Wigner function. In the linear limit, a quantum 
version of the Bohm-Gross dispersion relation was derived. Also certain nonlinear traveling 
wave solutions were obtained. 

It is the central purpose of this work, to extend the results of Ref. J6l to the 



electromagnetic case. Hence we define the moments 

n = Jdyf, (10) 

nu = Jdxfv, (11) 

Pij = m (^J d\ fviVj-nutu^j , (12) 

Qijk = m Jd\(vi-Ui)(vj-Uj)(v k -u k )f, (13) 

Rijki = m dv(vi-Ui)(vj-Uj)(v k -u k )(vi-ui)f, (14) 



and so on, as if / were a classical distribution function. Since all quantities are postulated 
in a gauge invariant way we can safely interpret n, u, Pjj etc. respectively as a particle 
density, a velocity field, a second rank stress tensor and so on. In particular, a scalar pressure 
p = (1/3) Pa and a heat flux vector g, = (1/2) Qjji can be deduced, where the summation 
convention is employed. Now the task is to obtain from the Wigner equation the equations of 
motion for the several moments, which will compose an infinite coupled hierarchy. 
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We also note that for the case of an isotropic distribution function, i.e. dependence of 
/ on the magnitude of the velocity only, all the odd moments must vanish from symmetry 
constraints, while the even moments are expressible in scalar quantities (by decomposition 
in terms of 5^). Moreover, for the case of local rotational symmetry, i.e. the existence of 
one preferred direction in (say, z) due to an external magnetic field or an initial temperature 
anisotropy, we have the the form 

P lJ =P ± h lJ + Ph l z J . (15) 

and 

Q lJ k = Q ± h {lk z k) + Qh l z J z k , (16) 

and similarly for higher order moments. Here we have introduced the projection tensor 
hij = 8ij —Ziij. These algebraic forms also solves the constraint equations (see below) that 
occur when assuming a stationary and homogeneous (but possibly anisotropic) equilibrium 
distribution^ 

For the sake of calculating the moments hierarchy equations, it is convenient to expand 
Av, B and E according to 

qh 2 eak , d 2 qh A Eak ^ d 4 

\2m i dvjdVm 540m 5 dVjdv m dv n dvi 

P _ E * 2 ,2 E d 2 h 4 4 d 4 

1 ~ 1 24m 2 J k 'dvjdvk 1920m 4 ^ kmn 1 d Vj dv k dv m dv n + ' 1 J 

h 2 d 2 h 4 d 4 

Bl = Bl ~ 2^ d i Bl dVJdV k + mjo^ 3 j kmnBi d Vj dv k dv m dv n + ' ' ' ' (19) 

disregarding higher order quantum corrections. The notation dj = d/drj is used whenever 
there is no risk of confusion. 

Assuming decaying boundary conditions, as far as the moment hierarchy is closed at the 
third-rank stress tensor, only the leading quantum corrections [the terms oc h 2 in Eqs. ([171 - 
(fT9l) l are needed. This is due to the structure of the higher order corrections. Indeed, these 
terms always involve at least fourth-order velocity derivatives and, for instance, 

d 4 f 

dyv i v i v k^-^ — 5-3- = 0. (20) 

dVaOVbdVcdd 

Therefore, only the semiclassical Wigner equation is needed, which does not mean that the 
quantum effects are necessarily small. It just happens that higher order quantum corrections 
would appear only for higher order moment evolution equations. 

Following Eq. ©, the semiclassical electromagnetic Wigner equation then reads 

^+v~ + i(E + vxB)~l /(r,v) 
at or m oy\ 

- d 2 E l I qh2£ijk d m B K I qh2£ 'J kVj d 2 B K d ' f (21) 

24m 3 i k 1 dvjdvjdvk 12m 3 m drjdvjdv m lArv? mn dvtdv m dv n 

$ We note that we can always decompose a moment of any order into its irreducible parts by picking an arbitrary 
direction and forming the projection operator orthogonal to that direction. 
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12m 4 



BidjBi 



B3B d "f ) 

dvidvjdvk 1 3 ' dvjdvkdvk J 



Notice that apparently the semiclassical electromagnetic Wigner equation, which has some 
interest in itself, was not discussed before in the literature. 
Calculating the moments, the result is 



Dn 

hnV-u 

Dt 



Dui 

~Dt 
Wj 

Dt 



DQi 



mn 



0. 



+ - (E + uxB). 

m 



— Pik^kUj Pj k d k Ui /jiV-U-|- £i mn Pj m B n -\- £i mn P[ m B n 

m m 

qh 2 qh 2 
+ yr-^m^i (ndjB k ) + j^EjkA - d k Q ijk , 



(22) 
(23) 

(24) 



Dt 



-Qijr u r U k — Qjkr d r Ui — Q kir d r U j — Qijk^ ' u — d r Rjj kr 



+ -!- (PijdrPkr + PjkdrPir + fly^r) + ~ (^irsQrjk + ZjrsQrki + ZkrsQrij) B s 

~ i€? ( d i Ek + d * Ei + d H + + 8 * 3i + 5kld ) Rl 



+ 



qh n 
12m 2 
qh 2 n 



(uxM-Kuxft + iux^B 



(25) 



Bin (djB r d s u k + d k B r d s u j) + £j rs (d k B r d s Uj + d}B r d s u k ) 

q 2 h 2 n 



12m 2 

T £krs 

(diB r d s Uj + djB r d s Ui) 



12m 3 



di(BjB k ) + dj(B k Bi)+d k (BiBj) 



When B = 0, Eqs. (1221 - 1251) recover the electrostatic equations [6J. In the limit h — > it 
reproduce the classical electromagnetic moment hierarchy equations [TT2l |T3l IT4TI . Quantum 
effects are explicit already in the transport equation for the pressure dyad, through the 
magnetic field. 

Previous approaches lfT5l derived quantum transport equations for charged particle 
systems assuming a local semiclassical Wigner function corresponding to a perturbed 
Maxwell-Boltzmann equilibrium. Here, however, the treatment include magnetic fields and 
is not semiclassical. A further approach for the derivation of quantum effects in macroscopic 
equations is through the eikonal decomposition of the wave functions of the quantum 
statistical ensemble and adequate simplifying assumptions ifToll . In both cases lfT5l [T6il the 
pressure dyad Pjj would be expressed as the sum of a classical part and a quantum part, the 
later one associated to a Bohm potential term in the force equation (|23l) . 

If we have used the gauge dependent Wigner function, it would not be possible to proceed 
exactly as in the classical case in the definition of the moments. Indeed, it would be natural to 
postulate them as 

~GD 



nu ■■ 



dp 



V-qA 



m 



f 



GD 



(26) 
(27) 
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q, = m Jdp & - -^ fM-nu^ , (28) 

@W = ^2 / d V(Pi ~ V Ai ~ mU ^> (P J ~^ A j~ mU j) (Pk - V A k ~ m Uk)f GV ■ (29) 

The same symbols n, u and Pij are used on purpose since Eqs. (|26l) - (|28l) produce the same 
expressions as from the GIWF, in spite of the fact that f GV itself is a gauge dependent object. 
However, from the equation satisfied by the usual Wigner equation |fTTfl one would obtain 



DPi i _ _ P 3 .. n 3 .. n T7 _ , 9„ n „ , 1 

s2 



Pik^kUj Pjk^kUi .P,jV • U -|- £i mn Pj m B n -\- £j mn Pi m B n 



D- t/v /v / in, »v t t / — 1 mwt /til — 1 

? m m 



-|^^ A - v "-^r, oo) 

containing gauge dependent quantum terms. The reason is that 

GST = G/y* - ^ (^5^^ + 4 A < + ^y) (3D 

is not gauge invariant. If Q9^ from Eq. (I3TI) is inserted into Eq. (|30l) one re-derive Eq. (T24l) 
for the pressure dyad on taking into account the Coulomb gauge which is assumed 0T1 in the 
evolution equation for f GD . 

Similarly the transport equations for the higher order moments are not gauge invariant. 
The conclusion is that to derive consistent equations from the usual Wigner function we would 
be obliged to modify the definition of moments. However, in this case there is the lost of one 
of the key advantages of using Wigner functions, namely the strict resemblance with the 
classical formalism. Also notice that if the heat flux tryad is set to zero the quantum term 
in Eq. (1301) is nonlinear for unmagnetized homogeneous equilibria, unlike Eq. (l24j) where a 
quantum contribution survives in this situation. 

In principle one could use the gauge dependent Wigner function to consistently calculate 
the higher order moments such as Qijk.Rijki an d so on. However, due to the fact that operators 
in quantum mechanics in general are non-commuting this cannot be done in practice. To see 
how this comes about we consider calculating the second order moment using the gauge 
dependent Wigner function. Calculating the second order moment Pjj(r,t) involves finding 
the expectation value of the operator, given bjQ 

\pi - qAi(r, t), [p j - qA j(t , t),8(t-r)] 



Uij ~4m 



(32) 



where [a, bj = db + bd denotes the anti-commutator. In order to calculate the expectation 
value using the Wigner formalism it is necessary to map the operator into a phase-space 
function using Weyl-correspondence ifTTll . This is done in practice by ordering the operators 
into a symmetric product of the position and momenta operators by using the commutation 
relations and then make the substitutions r — ► r and p — ► p. It turns out that the correct phase- 
space function is obtained by just making the substitution in the operator above without first 



§ The definition of the pressure operator in quantum mechanics is motivated by considering the Heisenberg 
evolution equation for the probability current operator which will be coupled to the divergence of the pressure 
operator. 
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Weyl order it. Hence we may calculate the pressure dyad using the gauge dependent Wigner 
function as 

Pij(r,t) = J [pi-qAtfj)] [pj-qAjir^t)] 8(r' -r)/ G V , V ,t) -mnu lUj (33) 

However, for the third order moment the correct phase-space function is not obtained 
simply by making the substitution r — > r and p — > p. Hence, calculating the correct third 
order moment using the gauge dependent Wigner function is complicated and involves Weyl 
ordering the corresponding operator so as to obtain the correct phase-space function. 

The GIWF has a modified Weyl ordering rule, discussed in fl4) and calculating the 
moments is done in complete analogy with the classical case, see Eqs. (fT0l)- (fT4)) . 



4. Transverse dispersion relation 

As an application of the fluid Eqs. ([22H25T ) we now consider linear transverse waves. 
Considering an one-component plasma, where the ions acts only as an homogeneous 
neutralizing background with number density no, the moment equations can be linearized 
around the equilibrium n = no, u = 0, Pjj = P-j , Qijk = 0, Riju = 0, E = 0, B = 0. To consider 
waves propagating in the z-direction with transverse polarization we let all fluctuations have 
the space-time dependence e lkz ~ imt and set E z = 0. Moreover we decompose the zeroth order 
pressure dyad as P-p = P±_(8j x 8j x + 8j y 8j y ) +P^8i z 8j z , where P± and are constants. 

It turns out that if we use the closure assumption R^ja = the quantum corrections to the 
transverse modes will not be retained so that to display the lowest order quantum corrections 
it is necessary to take into account also the contribution from the fourth order moment. As a 
closure assumption we use 

qh 2 / (Q) 3 (Q) 3 (Q) 3 (Q) 3 \ 
R m - ~ 4^3~2 [ F im d jkl + F jm d k\i + F km d lij + F lm d ijk ) E m , (34) 

adapted to the transverse wave case. The closure (l34l) is deduced systematically from the 
linearized equations satisfied by the fourth and fifth order moments, see Appendix A. Note 
that in principle the fourth order moment can have a nonzero equilibrium contribution 

^fjll ~ v t> wnere v t = ^J(2P±_ +P||)/(m«o) is the thermal velocity, but we will neglect this 
since we are looking only for the lowest order correction. Likewise for the terms ~ h 4 . Finally, 
it is worth to remark that in the classical limit the fourth order moment could be set to zero. 

The linearized equations can then be solved by first writing the magnetic field in terms 
of the electric field and then eliminating all quantities except the velocity so that we obtain 
the velocity in terms of the electric field. Coupling the resulting equation with Faraday's law 
via the current density J = qnoxx the dispersion relation 



ft) 2 - k 2 c 2 = (D 2 



k 2 P ± h 2 k 6 Pi 

1+ + - 



nornd) 2 4«om 3 ft) 4 



(35) 

is obtained. Here CO p = ^/noq 2 /(m£o) is the plasma frequency. If, instead, the closure 
Rijkl = was used, the term proportional to h would be absent in the dispersion relation. 
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In the simultaneous long wavelength and semiclassical limits, Eq. (1331) can be shown to 
admit an approximate solution 

+ : 5 \ - (36) 



, Pj.JP /j 2 fc 6 Pj_ 



To check the consistency, we need to compare to the results from kinetic theory. Here 
we are not concerned with Landau damping issues so that all integrals can be interpreted in 
the principal value sense. Assume 

E = Eiexp[/(Jfcz-a>f)]> (37) 

B = Bi exp[i(Jfcz - cat)] , (38) 

/ = /o(v) +/i(v) exp[i(*z - cot)] , (39) 

where k • E = as before and with the subscript 1 denoting first order quantities. The 
equilibrium Wigner function satisfy 



J d\f = n , J d\\f = 0. 



(40) 



Further we assume an equilibrium Wigner function such that fo = /o(vj_,v z ), where v\ = 
v 2 + v 2 . Notice that since there is no zeroth order magnetic field the perturbation velocity 
Av is also of first order. Hence Av does not contribute in the linearized Wigner equation ©. 
Using Eqs. dHHU) we get 



E = EL, B = BL, 

defining the operator 

sinh hk d 



L 



6 



2m <9i 



We note that 



dfo 



where k 



d V; 

kz. Moreover L 



m 

hk 
l 



/*k\ / hk 
/o v+— -f v- — 
V 2m \ 2m 



1 in the classical limit, since 



(hk d 
£ (2j + l)\ \ 2m~lh z 



2j 



l fhk\ 2 d 7 



l + ^T — 



24 \ m J dv 



+ ■ 



(41) 



(42) 



(43) 



(44) 



Then linearizing the Wigner equation © and from the Maxwell equations with charge 
and current densities q (f d\f — no) and q f dwf respectively, the result is 

v] Lfo 



k 2 CO 2 
co 2 = co 2 + c 2 k 2 + - 

1 2n 



V 



d\- 



(45) 



(co-k-v) 2 ' 

where c is the speed of light and co p is the plasma frequency. In comparison to the classical 
transverse dispersion relation, the only change is the replacement fo — > fo = Lfo- In a 
classical picture it is as if the particle velocities were reorganized through the diffusive 
operator L. Also notice that still fo = /o( v _i_5 v z) ■ Moreover, the quantum diffusion induced by 
the operator L preserves the number of particles, since / d\fo = f d\ fo due to Eq. (|44|) 
under decaying boundary conditions. Figure 1 shows the effect of L on the equilibrium 
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v 



Figure 1. Quantum diffusion on the equilibrium Wigner function fo = 

/r(vj_) exp[— v^/(2vq)]. Here f\\ = L(exp[— v^/(2vq)]). Values of the parameter 
H = hk/(2mv ) are// = 0,1 and 2, so that /j| (0) = 1, 0.86 and 0.60 respectively. 



/o = /r(vj_) exp[— v^/(2vq)], for different values of the non dimensional parameter H — 
hk/ (2mvo). In the simultaneous long wavelength and semiclassical limits and retaining only 
the leading ~ v\ thermal corrections, Eqs. (1361 ) and (1431 ) give the same result via the natural 
identification P± = (m/2) J dw^fo. This conclude the equivalence between the moments and 
kinetic theories, in the fluid limit. 

To compare, the transverse dispersion relation following from the gauge dependent 
Wigner function [HUES can be expressed as 



CO 



C0p+c 2 k 2 



mCOp 
2noh 



dv- 



co-kv 



fo 



-fo 

2m 



hk 

2m 



or, using Eq. (1431) . as 



2 1 « 2 k 2 



CO 



2 k 



co p + c 



2 n Q 



d\- 



dfo 



,(46) 



(47) 



CO — k ■ v \dv z 

An integration by parts then shows the equivalence with the gauge invariant transverse 
dispersion relation Eq. (1431) . Therefore the gauge choice issues tend to be crucial only for 
the nonlinear regimes, as also manifest in the gauge dependent nonlinear term in Eq. (1301 for 
the pressure dyad. However, in the case of non-homogeneous equilibria the use of a gauge 
independent electromagnetic Wigner equation is advisable even for linear waves. 



5. Conclusion 



The moment hierarchy equations derived from the GIWF electromagnetic evolution equation 
is obtained. The advantages over the gauge dependent Wigner formalism are stressed. 
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Discrepancies tend to be prominent in the nonlinear regimes and for higher order moments of 
the Wigner function. The fluid-like equations (I22l)-(l25l). closed at the transport equation for 
the heat flux tryad, is applied to the propagation of linear transverse waves. Good agreement 
is found when comparing with the results from kinetic theory, in the long wavelength 
approximation. A key ingredient to a successful macroscopic theory is an adequate closure 
of the moment equations and a recipe for solving this question is proposed, see Appendix A. 
The approach is not restricted to particular local equilibrium GIWFs and is not based on a 
Madelung decomposition of the quantum statistical ensemble wave functions. The moment 
equations (l22~l)-(l25T) is an adequate starting point for studying the nonlinear aspects of quantum 
plasma problems involving magnetic fields, e.g. via numerical simulations. 
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Appendix A. The closure problem 

The closure (|34l) can be deduced systematically from linearized higher order moment 
equations. Let Sukim De the fifth order moment defined in analogy to the third and fourth 
order moments, see Eqs. (fT3l4T4l) . The sixth order moment will be set to zero. The evolution 
equations for the fourth and fifth order moments are derived following the same steps as when 
Eqs. (f2~2l - [2~5l) were derived starting from Eq. ([2T]) . Since they are quite complicated we here 
only include the linear terms, which gives 




12m 3 [ mm \ J" d kl+r kn o-ij+r jn dj k +r jk d ln + f u O jn + ^j d kr 

+ t -jnmykn a li^Un a ik + ^in "kl^^kl a in + Ui a kn^ Uk "in) 

4_ p, ( p(«) ? 2 + p(0) -a , p(0) -a p (o) -a p (o) ?2 p (o) ?2 \ 

+ £lnm{t J in d jk + F jn d ki + F kn d ij +F ij d kn+ F jk d in+ F ki 6 }n)\ H m 



) 




(A.1) 



(A.2) 



Moment hierarchy from quantum kinetic theory 



12 



After Fourier transforming and inserting Syjfc/m from Eq. (IA.2I ) into Eq. (IA.1I) , Eq. 
(|34|) is derived using Faraday's law. The procedure is adapted to the present equilibrium 
(homogeneous, no streaming particles, no heat flux, negligible higher order thermal effects). 

It turns out that due to a cancellation arising from Faraday's law in the transverse case, 
the result in Eq. (|34|) is correct even if Eqs. (IA.H - IA.2l) were extended to include ~ h 4 
terms. To obtain the next order quantum effects dispersion relation using the fluid theory, it 
is hence necessary to include higher order moments. In this example, the sixth order moment 
disregarded in Eq. (IA.2I) . 

From the above we can infer a general recipe for the closure of the fluid-like system up to 
the N th — moment: Fourier transform the linearized evolution equations for the (N+l) th and 
(N + 2) th moments, setting the (N + 3) th moment to zero. In this way we derive an expression 
for the (N + l) th moment, so as to close the system for the N moments. The form of the 
linearized equations depends on the particular equilibrium. Naive closures like setting the 
(N + \) th — moment directly to zero tend to produce fake results when comparing to kinetic 
theory. This is in sharp contrast to the simplicity of the electrostatic case, where faithful 
equations are obtained already defining the fourth order moment to be zero (H. 
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